Astronomy & Astrophysics manuscript no. 


February 2, 2008 


(DOI: will be inserted by hand later) 





Multi-wavelength study of the G 82.2+5.3 supernova remnant 

F. Mavromatakis^ , B. Aschenbach^, P. Boumis^''^, and J. Papamastorakis^'^ 



University of Crete, Physics Department, P.O. Box 2208, 710 03 Heraklion, Crete, Greece 
. . Max-PIanck-Institut fr extraterrestrische Physik, Postfach 1312, D-85741 Garching, Germany 

CO . Foundation for Research and Technology-Hellas, P.O. Box 1527, 711 10 Heraklion, Crete, Greece 

o ■ 

^ i Received 8 August 2003/Accepted 6 November 2003 

^ I Abstract. We present the first CCD flux-calibrated images of the supernova remnant G 82.2+5.3 in major optical 

O ■ emission lines. The medium ionization line of [Olll]5007 A provides the first direct evidence of optical emission 

' originating from G 82.2+5.3. Filamentary emission is detected in the west and east areas of the remnant, roughly 

I , defining an ellipsoidal shell. The [O ill] emission is rather well correlated with the radio emission suggesting their 

■ association, while typical fluxes are found in the range of 20-30x10"^^ erg s~^ cm~^ arcsec"^. Deep long-slit 

spectra taken at speciflc positions of the remnant verify that the detected filamentary emission originates from 
, shock heated gas, while the diffuse [O ill] emission in the south results from photoionization processes. The spectra 

^ ' further suggest shock velocities around 100 km s~^ and low electron densities. The X-ray surface brightness is 

I quite patchy, missing obvious limb brightening and is dominated by a bright bar-like emission region which is 

CO . off-set from the geometric center by ~9'. The X-ray emission is thermal and requires two temperatures of 0.2 

^ • keV and 0.63 keV. The bright bar region shows overabundant Mg, Si and Fe, which might indicate still radiating 

ejecta matter. The azimuthally averaged radial surface profile is consistent with the matter density changing 
with distance r from the center oce~'^'° with a characteristic angular length of 36', or, alternatively, with an 
r~^/^ density profile. The matter inside the remnant is quite likely structured like a porous cloudy medium. The 
average matter density is ~0.04xd^g'^ with di.e the distance in units of 1.6 kpc. Because of the low density and 
• the long cooling times involved the remnant is more likely to be in the adiabatic phase, which is consistent with 

the densities derived for the X-ray plasma and the optical line emission, but it is not excluded that is has reached 
Q . the radiating phase. This, however, would imply a lower density, greater age and much larger distance, at the 

' edge of the upper limits obtained from Nh and E(B-V). 
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1. Introduction ages and the radio maps of G 82.2+5.3. Optical spectra 

■ taken from Sabbadin 1)1976(1 show the characteristic sig- 

, The supernova remnant G 82.2+5.3 (W 63) is found in ^^^^^^ of emission from shock heated gas ([S ii]/Ha^0.5). 

the complex region of the Cygnus constellation. Radio j^^^^^^ ^ Gonzalez mD estimated expansion velocities 

observations at several frequency bands established the ^^^^ ^p^^^^l filaments of ^35 km s"! and -70 

nature of the object through its non-thermal radio spec- ^-i ^ depending upon the size assumed for the remnant. 

trum (Velusamy & Kundu 19741 Angerhofer et al. 119771 rn, v j j- i i • i 

^ / ■ — ' ° ' ' ihe X-ray and radio morphologies oi the remnant are 

Wendker Ml^- The remnant displays a non-circular, ^^.^^ different. Rho & Petre 099^ proposed that it be- 

elliptical shape of —70 xlOO in the radio map at 11 , ; .i i r • j u i ^ t j j 

, . n 1 • 1 longs to the class ot mixed-morphology remnants. Indeed, 

cm. The distance to the remnant is not well determined ,i v ... .,, , , . ,., 

the A-ray emission is neither shell-like nor plerion-like 

and distances in the range of 1.3-1.9 kpc have been re- , , , • r ii v • ■ ^ r 

f , . ^ but the vast majority oi the A-rays originate trom areas 

ported (Rosado & Gonzalez 119811 and references therein). juj-j. ij 4- ^ j^i, 

f „ ^ „, , ~ — T „ , , ■ , around but not peaked on the remnant center, and there 

Interference-filter photographs from Parker et al. 11197911 . , . r u • uj- • r^i, ^- v 

, , 1 rn 1 1 — IS obvious limb brightening, ihe entire A-ray region 

revealed bright extended structures in Ha and S II . There . t,i i, j u 4-1, j- t, n oi, c n j- 1 ^ AA6 I1 

° ' IS roughly bound by the radio shell. Rho & Petre 1)199811 

is no obvious correlation between the low ionization im- , . -Dr\a\rr jj-j -uj-i, j-- ^ i, 

analysing ROSAi data describe the entire spectrum by a 

Send offprint requests to: F. Mavromatakis,e-mail: thermal, single temperature spectrum with kT -0.2 keV 

f otisOphysics uoc gr Subjected to interstellar absorption with a column density 

* Present address:Institute of Astronomy & Astrophysics, of -4 x lO^i cm-^ although with large uncertainties. 

National Observatory of Athens, I. Metaxa & V. Pavlou, P. Flux calibrated images in major optical emission lines 

Penteli, 15236 Athens, Greece were obtained and deep long-slit spectroscopy of selected 
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Table 1. Journal of the Observations 





Ac 


AA 




Exp. times'' 


Filter 


(A) 


(A) 


Date (UT) 


(No of frames'") 


Ha + [N II] 


6560 


75 


25-26/07/2001 


7200 (3) 


[Sii] 


6708 


20 


25/07/2001 


7200 (3) 


[Oiii] 


5005 


28 


26/07/2001 


9600 (4) 


[On] 


3727 


25 


26/07/2001 


7200 (3) 


Cont red 


6096 


134 


25-26/07/2001 


720 (8) 


Cont blue 


5470 


230 


26/07/2001 


360 (2) 



Total exposure times in sec 
^ Number of individual frames 



areas of interest was performed for a detailed study of the 
optical emission in the area of G 82.2+5.3. In addition, we 
analyzed ROSAT All-Sky survey and re-analyzed ROSAT 
pointed data as well as ASCA GIS and SIS data from one 
pointing. Information about the observations and the data 
reduction is given in Sect. 2. In Sect. 3 and 4 we present the 
results of our imaging and spectral observations. Finally, 
in Sect. 5 we discuss the properties of the remnant and its 
environment. 

2. Observations 

2.1. Optical images 

The current wide field observations were performed with 
the 0.3 m Schmidt-Cassegrain telescope at Skinakas 
Observatory, Crete, Greece. The complex field of 
G 82.2+5.3 was observed on July 25 and 26, 2001. The 
1024 X 1024 SITe CCD used during these imaging ob- 
servations resulted in a 89' x 89' field of view and an 
image scale of 5" per pixel. A journal of the observations 
together with information about the filters are given in 
Tabled The images presented in this work are the aver- 
age of the available frames from the individual filters. The 
astrometric solutions for all data frames utilized the HST 
Guide star catalogue (Lasker et al. 11999(1 and all coordi- 
nates quoted in this work refer to epoch 2000. All frames 
were projected to a common origin on the sky. 

Standard IRAF and MIDAS routines were employed 
for the reduction of the data. All frames were bias sub- 
tracted and fiat-field corrected using a series of well ex- 
posed twilight fiat-fields. The spectrophotometric stan- 
dard stars HR5501, HR7596, HR7950 and HR9087 were 
observed in order to obtain absolute flux calibration 
(Hamuy et al.^M^^- 

2.2. Optical spectra 

The 1.3 m Ritchey-Cretien telescope at Skinakas 
Observatory was used to obtain long-slit spectra on June 
21, 22, and July 21, 22 2001. The data were taken with a 
1300 line mm~i grating and a 800 x 2000 SITe CCD cov- 



Table 2. Spectral log 



Slit centers 




Exp. times* 


a 


S 


(No of spectra*") 


20"21™39'.3 


45°59'30".l 


16500 (5) 


20'' 20^04" .0 


44°54'48".5 


7200 (2) 


20" 16™ 15". 6 


45°51'58".9 


7200 (2) 


20" 16"' 15' .4 


45°58'01".4 


7200 (2) 



* Number of spectra obtained 
*" Total exposure times in sec 



ering the range of 4750 A - 6815 A. The slit had a width 
of 7'.' 7 and a length of 7'. 9 and, in all cases, was oriented in 
the south-north direction. The coordinates of the slit cen- 
ters, the number of available spectra from each location 
and the total exposure time of each spectrum are given in 
Table El The spectrophotometric standard stars HR5501, 
HR7596, HR9087, HR718 and HR7950 were observed for 
the absolute calibration of the spectra of G 82.2+5.3. 

2.3. The ROSAT and ASCA observations 

The Cygnus region was scanned during the ROSAT 
All-sky survey in early December 1992. X-rays from G 
82.2+5.3, which is located at the northern tip of the so- 
called Cygnus Superbubble shell (Uyaniker et al. I2001|l 
were clearly detected and an image of the remnant and 
the environment could be constructed. The net observ- 
ing time was ^^1200 s. The remnant under study was 
also observed several times with the ROSAT PSPC in 
pointing mode, the longest observation of which was 
for ^^8 ks on 26 October 1993 aiming at the center of 
the remnant (sequence #500286p). Additional ROSAT 
observations were targeted at the west of the remnant 
(#500219p/500219pl), the north-east (#500217p) and 
the south-east (#500218p). Apart from the pointing to- 
wards the centre only the south-east pointing provided 
sufficient counts to create a spectrum over the central 40 
arcmin of the PSPC. We also found one observation in the 
ASCA archive and both the GIS and SIS data could be 
used for a spectral analysis. 

3. The optical emission 

3.1. The Ha + [N ii] and [S ii] line images 

The images in Ha + [Nil] (Fig. P) and [Sii] (Fig. ^ of 
the remnant look complex due to the presence of sev- 
eral bright and dark nebulae. The dark nebulae LDN 908 
and 910 (Lvnds [T962|l are found in the south areas of our 
field, while LBN 081.88+04.79 (LBN 314; LyndsTMS') is 
located in the south-east and LBN 082.79+05.83 (LBN 
325), LBN 083.03+05.78 (LBN326) located to the north 
may contribute to the emission in our field. Lynds H1965|l 
also reports LBN 081.62+05.47 (LBN 310) as a very ex- 
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Fig. 1. The field of G 82.2+5.3 in the Ha + [Nil] filter. The image has been smoothed to suppress the residuals 
from the imperfect continuum subtraction. The open and the filled circles designate the positions of bright and dark 
nebulae, respectively, while the numbers in parentheses represent the surface of each nebula in square degrees as given 
by Lynds H1962I [19651. The open polygons show the positions where emission line ratios are measured (see also §3.3). 
The shadings run linearly from to 400 x 10~^^ erg s~^ cm~^ arcsec"^. The line segments seen near over-exposed 
stars in this figure and the next figures are due to the blooming effect. 



tended Hll region which seems to overlap our field. The 
complexity of the environment of the remnant does not al- 
low the direct identification of any optical emission match- 
ing the morphology of the known radio emission. However, 
since our images are flux calibrated they can be used to 
approximately determine the emission characteristics. The 
ratio of the low ionization images of [S ii] and Ha + [N ii] 



shows that the bulk of the emission in the west results 
from photoionization, while the weak emission seen to the 
south of the bright star GSC 3572-2102 probably orig- 
inates from shock heated gas. In this area, roughly de- 
fined along a ~ 20^22^^30" and between S ~ 45°47' and 
45°12', we estimate a [Sii]/Ha ratio of ~0.7-0.8, while 
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Fig. 2. The [S ii] image of the area around G 82.2+5.3. The morphology in these emission Unes is similar to that of the 
Ha + [N ii] lines. The image has been smoothed to suppress the residuals from the imperfect continuum subtraction, 
while the shadings run linearly from to 60 x 10~^^ erg s~^ cm~^ arcsec"^. 



LBN 082.79+05.83 and the bright structure in the north- 
west are characterized by ratios very close to ~0.4. 

3.2. The [0 ii] and [0 iii] images 

The [On] image of the area around G 82.2+5.3 is not 
shown here since its morphology is generally similar to 
that of the Ha + [Nil] and [Sii] low ionization images. 
Emission from the Hii region in the [Oii]3727 A line is 
somewhat suppressed but still the shape of the remnant is 



not clearly outlined. No immediate correlation of the op- 
tical emission with the known radio emission can be seen. 
In strong contrast to the low ionization images, the field 
of G 82.2+5.3 appears markedly different in the medium 
ionization line of [O ill] 5007 A (Fig.|21). The overall emis- 
sion is significantly attenuated and the image appears to 
be free of complex structures. Optical radiation possibly 
from LBN 082.79+05.83 and/or 083.03+05.78 is seen in 
the north edge of the field of view. The [O ill] image reveals 
filamentary emission in the west and east areas as well as 
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Table 3. Relative line fluxes 





I (west) 


II (east) 


III (south) 


Line (A) 








4861 H/3 


25 (19) 


15 (10) 


11.4 (66) 


4959 [Oiii]i 


34 (20) 


159 (104) 


1.5 (11) 


5007 [Oiii]2 


105 (57) 


511 (293) 


5.0 (36) 


5872 He I 


— 




3.5 (37) 


6300 [Oi]i 




75'' (32) 




6364 [O1J2 




22 (IS) 




r /I Q fNT ttI 
D04o [IN IIJl 


io [i-o) 


A A i' A A\ 

44 (44 ) 




ooDo xlQ: 




inn ^'7'i\ 
lUU [10) 


inn ('c;on^ 
lUU [0Z[)) 


«KiiA nS,T ttI 

D0o4 [i\ nj2 


Do (^00 ) 


1/11 ^^ ^ ^\ 
i4i t^iioj 


61. L yZKyO) 


DDio xlei 






1 Q /I ''7\ 

1.3 (17) 


fi7i (s Till 




i y i J 


ins ^^^'\'\ 


6731 [Sii]2 


25 (29) 


56 (53) 


7.6 (82) 


/\Dsoiute tict nux 


1 n 


4.0 




TT „, /TT /□ 

net/ tip 


4.U l^ioj 




0.0 I^DOj 


[u 111]/ np 


O.D l^J-oJ 


('1 n^ 




[Sii]/Hq 


0.6 (39) 


1.34 (56) 


0.18 (44) 


I(6716)/I(6731) 


1.40 (23) 


1.39 (43) 


1.42 (66) 




0.38 [0.07] 


1.1 [0.1] 


1.41 [0.01] 


E(B-V) 


0.25(0.05] 


0.73[0.07] 


0.94[0.01] 


He 


< 80 cm"'' 


< 30 cm"^ 


< 30 cm"^ 




< 100 km s"^ 


> 100 km s"^ 





^ Uncorrected for interstellar extinction 

^ Listed fluxes are a signal to noise weighted average of the individual fluxes 

^ Present in only one of the spectra 

* In units of 10~^^ erg s~^ cm~^ arcsec"^ 

^ The logarithmic extinction is determined as c = 1/0.331 • log((Ha/H/3)obs/3) 
Numbers in parentheses represent the signal to noise ratio of the quoted fluxes 
The la error is given in square brackets for the logarithmic extinction 
All fluxes normalized to F(Ha)=100 

Table 4. X-ray spectral fits 



Instrument 


Nh 

lO^icm-2 


kT 
keV 


EM 

lO'^cm"'^ 


Xrcd 


dof 


GIS/SIS 

PSPC^ 

PSPC^ 


^- ' -0.4 

4.0+1:° 


0.63+-.- 


1 5Q+0.20 

0.57i°:^° 
4.70^^:1 


0.99 
0.74 
0.79 


198 
22 

48 



^: one temperature vmekal model 

^: two temperature vmekal model; values of low temperature 
component listed, values of high temperature component 
restricted to GIS/SIS fit range. 



diffuse emission. Bright filaments are detected along a ~ 
20'^16'" and 20^22°" and around 6 ~ 45°40'. These two fil- 
amentary structures define rather well the opposite sides 
of an ellipsoidal shell which is in agreement with the mor- 
phology of the radio emission. In fact the correlation of 
the [O III] emission in the east with the radio emission ob- 
served at 326 MHz and 4850 MHz (Fig. gj) suggests its 
association to G 82.2+5.3. Fainter structures are also de- 
tected in the south and inner regions of the field. Faint 



emission from a region shaped as an arc is detected to the 
east of GSC 3572-2102, roughly along a ~ 20^22'"40^ 
from S ~ 45°30' to 46°30'. This faint filamentary struc- 
ture is not continuous but gaps in intensity are detected 
along its extent. This may indicate the presence of inho- 
mogeneities in the interstellar "clouds" resulting in strong 
variations of the shock velocity. It is known that the shock 
velocity greatly influences the [O ill] flux (Cox & Raymond 
11985(1 . The typical projected angular width of the [Olll] 
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filaments is ~30", equivalent to 0.2 pc at an assumed av- 
erage distance of 1.6 kpc (see §1). 

3.3. The long-slit spectra from G 82.2+5.3 

The deep low resolution spectra were taken in the vicinity 
of the bright stars GSC 3572-2102 and GSC 3572-2103, 
as well as in the south in an area where strong radio emis- 
sion overlaps with strong Ha -I- [N ii] and diffuse [O ill] 
emission (Fig.^). The optical spectra from areas I and II 
(Fig. 01 indicate emission from shock heated gas ([S ii]/Ha 
> 0.4), while that from area III seems to originate from 
photoionized gas (Table This is because we cannot ex- 
clude the possibility that the spectrum from area III may 
be contaminated by emission from a nearby Hii region, 
thus causing the low [Sii]/Ha ratio of 0.18. 

Due to the presence of strong Hii emission around 
area I, we report the procedure followed with the spec- 
tral analysis. The first set of spectra were taken with the 
slit crossing the [Olll] filament. However, the area along 
the slit and in the north of the [O ill] filament contained 
very strong H ii emission. In order to take into account the 
diffuse emission present in the field, the second set of spec- 
tra were taken with the slit shifted to the north by 6'. In 
this way the area around 5 ~ 46°00', containing relatively 
weak Hii emission, was within the slit and was used as 
background. It is clear that the spectral results from area 
I should be used under the assumption of constant back- 
ground intensity. The low [O ill] /H/3 ratio of ^^6 implies a 
shock speed ^^100 km s~^ or less, while the sulfur line ra- 
tio suggests a pre-shock cloud density of a few atoms per 
cm^ (e.g. Raymond et al. EHHl Osterbrock UnHHIl ■ The 
spectrum from area II suggests electron densities below 
~80 cm~^, while the strong sulfur line emission may sug- 
gest a largely neutral pre-shock medium. In addition, this 
spectrum displays very strong [O ill] emission relative to 
H/3, a characteristic of shocks with incomplete recombi- 
nation zones (Raymond et al. ll988|l . This high [O III] /H/3 
ratio of ~45 further suggests shock velocities greater than 
~100 km s-^ 

The Ha/H/3 ratio of 4.0 (±0.2) from area I results in 
a logarithmic extinction of 0.38 (±0.07) or an E(B-V) of 
0.25 (±0.05), while in area II the logarithmic extinction of 
1.1 (±0.1) is equivalent to an E(B-V) of 0.73 (±0.07). The 
logarithmic extinction is calculated through the relation 
c = 1/0.331 • log((Ha/H/3)obs/3) and uses the interstellar 
extinction curve of Kaler H197f)|l as implemented in the 
nebular package (Shaw and Dufour'1995') within the IRAF 
software. The E(B-V) is calculated with the aid of the 
relation E(B-V) = 0.664 ■ c (Kaler [W7^ Aller lTOMjl . We 
note that the signal to noise ratios quoted in Tabled do 
not include calibration errors which are ~7%. 

4. X-ray emission from G 82.2+5.3 

The field of G 82.2+5.3 was observed by ROSAT in 
pointed mode on October 26, 1993 (ROR 500286p) for 
8 ks. The image of this observation is shown in Fig. El 



On purpose we show just the count/pixel image and did 
not correct for exposure, to illustrate where the remnant 
is highly underexposed due to the artifacts of the PSPC, 
clearly seen as a dark bar running from the lower left to 
the upper right of the image and other ribs, as well as 
a dark ring 40 arcmin wide around the center. Since the 
image is not exposure corrected the intensity close to the 
edge is brighter by a factor of about two than the bright- 
ness in the image indicates. The brightness appears quite 
clumpy and is dominated by a relatively bright bar off- 
set from the geometric centre to the north by about 9 
arcmin. The bar measures about 24 arcmin in the north- 
south direction and about 18 arcmin along the east-west 
direction. In Fig. the PSPC image is shown along with 
contours of the radio emission at 326 MHz, which best 
define the elliptical shape of G 82.2+5.3. It is also evident 
from this figure that the X-ray emission from the rem- 
nant does not fill the PSPC field of view, which is sup- 
ported by the all-sky image on an even wider scale. Thus 
it is possible to determine the background from the same 
observation. Checks with the background taken from the 
adjacent pointings have been made and the source spec- 
trum does not change. G 82.2+5.3 was also observed by 
the ASCA satellite on June 12, 1997. The ASCA X-ray 
telescope was pointed to the central area of G 82.2+5.3 
for a total of 82 ks. Given the smaller field of view of the 
CIS and SIS detectors compared to the ROSAT PSPC, 
this data can only be used for spectral analysis (see §4.2). 

4.1. Image analysis 

The first and so far only results of an analysis of the 
ROSAT data of G 82.2+5.3 were reported by Rho & Petre 
|1998), but no images were shown. The X-ray image looks 
like an SNR broken into pieces resulting into a patchy 
appearance. The central bar is the brightest component 
of the remnant and has led Rho & Petre 0i398) to put 
G 82.2+5.3 in the class of mixed-morphology remnants. 
First, we have searched the remnant for point sources. 
The ROSAT source detection algorithm reveals the exis- 
tence of 31 point sources with a likelihood > 10 or 4.1 a 
across the face of the remnant. Within a circle of 40' di- 
ameter centered on the apparent SNR centre 17 point-like 
sources are found out of which four sources do not have 
optical counterparts down to me < 22, with separations 
from known optical sources of >10". The source existence 
significance of 4.1 a corresponds to a source count rate of 
0.01 counts s~^ which corrected for an absorption column 
density of Nh = 4.7x10^^ cm~^ (§4-2) corresponds to an 
unabsorbed (0.2 - 2.4 keV) energy flux of ~ 2xl0~^^ erg 
cm~^ s~^ for a blackbody source spectrum of kTbb = 100 
eV. Of course, with a detection significance of 4.1 a it can 
not be excluded that these point-like sources are actually 
statistical fluctuations above the diffuse remnant emission. 

Although the surface brightness is quite low, we can 
integrate over the azimuth and construct its radial pro- 
file including every PSPC count independent of energy. 
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Fig. 3. The medium ionization line of [Oiii]5007 A provides the sharpest view of G 82.2+5.3 (left). The image has 
been smoothed to suppress the residuals from the imperfect continuum subtraction and the shadings run linearly from 
to 30 X 10~^^ erg s^^ cm^^ arcsec^^. The long rectangles indicate the projected positions of the slits on the sky 
and the numbers designate the areas discussed in the text. The right figure shows single long-slit spectra from areas 
II and III. 




Fig. 4. The correlation between the [O ill] emission and the radio emission at 326 MHz (left) and 4850 MHz (right) 
is shown in this figure. The 326 MHz radio contours (Rengelink et al. I19H7|I scale linearly from 0.01 Jy/beam to 0.09 
Jy/beam, while the 4850 MHz contours scale from 5 10^^ Jy/beam to 0.2 Jy/beam (Condon et al. 11994(1 . The optical 
and radio data are rather well correlated in the east areas of the remnant. The radio emission in the west seems to be 
weeker and mainly related to the H ii region. 
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Background was subtracted, the profile was corrected for 
vignetting and tlie integration was done over 360° rather 
than just 40° which was used by Rho & Petre (11)98). 
Furthermore we centered the profile on the apparent geo- 
metric centre (a ~ 20'"19™07% S ~ 45°31'12") of the rem- 
nant rather than on the brightest patch (c.f. O. It turns 
out that the density profile can be well described by an 
exponential law (oc e""'^''") with a characteristic angular 
length of 36'(±3') up to a distance of 45'.5 (±1'.5). At 
larger distances the density can be described by the same 
exponential law but enhanced by a factor of 1.41q 2, which 
might be interpreted as an indication of limb brightening. 
The brightness profile is also consistent with a density 
power law, such that the density scales inversely with the 
square root of the distance from the center. The actual 
data allow us to trace the density profile out to an an- 
gular distance of ~50', which we assume to be the radius 
of the remnant (Fig. EJ. In summary, the average matter 
density drops by just about a factor of four from the cen- 
ter to the edge. Of course, the radial brightness change 
could also partially be produced by a radially changing 
filling factor. The patchiness of the surface brightness dis- 
tribution is then a matter of local density variations and 
volume filling factors. This result is quite in contrast to 
the result of Rho & Petre EMI who, however, have cho- 
sen a fairly small sector of the remnant (1/9 of the whole 
perimeter) and probably have centered the profile on the 
brightest spot. 

4.2. Spectral analysis 

In order to study in as much detail as possible the spectral 
properties of G 82.2-1-5.3, we analyzed archival ASCA data 
and the ROSAT data. Due to the low surface brightness 
of the object, both the ASCA and ROSAT PSPC data are 
characterized by low signal to noise, not allowing a spa- 
tially resolved spectral analysis. Thus the analysis of the 
GIS data was restricted to the brightest region within a 
circle centered on a ~ 20^19'"14% S ~ 45°47'00" and 7'.4 
radius. However, the center of the ASCA SIS detectors is 
offset by ~6' with respect to the center of the GIS detec- 
tors. Since the SIS data partially overlap the GIS data, 
counts from a 3' x 3' square area were extracted from the 
former detectors. 

A non-thermal component is immediately ruled out 
by the ASCA data since a power-law fit is not accept- 
able (reduced x^=1.74 for 197 degrees of freedom). A 
mekal model fitted simultaneously to the GIS/SIS data 
with non-cosmic abundances results in a reduced of 
1.15 for 196 degrees of freedom (dofj, which is formally 
acceptable. However, the presence of a systematic trend 
in the residuals in the 0.9 to 1.9 keV range suggests its re- 
jection, while including a power-law component does not 
lead to an improved fit. The GIS/SIS data were also fit- 
ted with a vmekal model and a reduced of 0.99 for 198 
dof (Fig. d) was obtained. The best fit parameters of this 
model are given in Tabled Adding a second component 



(power-law) does not substantially improve the fit. The 
vmekal fits also constrain the abundance of O (1±1), Ne 
(0.lt°:I), Mg (2.8lJ:37), Si (4.0^1:^) andPe (3.512J); num- 
bers in parentheses are abundances relative to solar values 
with la errors, which shows that Mg, Si and Fe definitely 
differ from solar values. The Si line is very clearly seen in 
both the GIS and SIS spectra (c.f. fig. Ej) at around 1.8 
keV. This is interesting because it might mean that the 
bright patch, i.e. the near center component, is dominated 
by ejecta matter rather than circumstellar or interstellar 
matter. For the ROSAT analysis we used the same field 
and the same abundances with a one temperature vmekal 
spectrum. Fig.|Slshows the error contours relating the tem- 
perature kT and the interstellar absorption column den- 
sity Nh. As the contour lines show there is basically no 
overlap between the ASCA and the ROSAT fits, except 
at the 4(7 level (99.99% probability) for each of the two 
instruments. Clearly a one temperature model does not 
represent the spectrum. Obviously a component with a 
temperature too low to be observed by ASCA - we have 
restricted the ASCA analysis to energies > 0.8 keV - has 
to be added which on the other hand is dominating the 
ROSAT spectrum, or vice versa. So we fitted the PSPC 
spectrum with a two component vmekal model, restricting 
the values of the high temperature, the associated emis- 
sion measure as well as Nh to the range of ±1ct around 
the best fit found in the ASCA data modeling. The best 
fit values for the low temperature component are given in 
Table0]as well. In summary the high temperature is found 
to be around 0.6 keV and the low temperature around 0.2 
keV. The emission measure of the low temperature com- 
ponent is about three times that of the high temperature 
component. Both the Nh and the value of the low tem- 
perature are consistent with the results of Rho & Petre 
119981 but the existence of a high temperature component, 
already indicated by a deeper analysis of the PSPC spec- 
trum, is very clearly demonstrated by adding the ASCA 
data to the analysis. 

The emission measures of both components allow us to 
estimate the corresponding densities of the X-ray emitting 
gas. It is found that, for a distance of 1.6 kpc, the aver- 
age hydrogen number density of the harder component is 
~0.03 cm~^, while that of the softer component is ^0.05 
cm^'^ and both scale inversely with the square root of the 
distance to the remnant. These numbers adopt an emission 
volume of 7r-(7.4 arcmin)^-100 arcmin, so that the hydro- 
gen number density in cm^'^ is uh ~ 0.657 x -\/ (EM/di.e), 
with EM the emission measure (c.f. Table^J). We note here 
that a change by a factor of 4 in distance would alter the 
densities by only a factor of 2. It is clear that the X-ray 
gas is quite thin and that the remnant expands in an even 
thinner medium, if we consider the matter inside as being 
compressed by shocks. 

As mentioned above another pointing on G 82.2-1-5.3 
with spectrally useful data has been performed. The spec- 
tral analysis of the PSPC pointing centered on the south- 
east of the remnant shows that the spectrum is fully con- 
sistent with that of the central part, and no significant 
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differences can be established. In fact the fit with a one 
temperature vmekal model produces the same parame- 
ter values within la error bars. The existence of a high 
temperature component cannot unambiguously be demon- 
strated, it might be present or might be missing. ASCA 
data on the field are not available. 

5. Discussion 

The supernova remnant G 82.2-f 5.3 is found in a complex 
area in the Cygnus constellation where strong Hll emis- 
sion dominates the red part of the optical spectra. Most 
of the radio emission is detected in a shell which basically 
defines the boundaries of the remnant. However, a defi- 
ciency of radio emission at the west part of the remnant 
is present, at least in certain radio frequencies as well as 
in the X-rays. In this work we present the first flux cali- 
brated CCD images of G 82.2-1-5.3 in the emission lines of 
Ha-I- [N ii] , [S ii] , [O ii] , and [O ill] . In addition, results from 
a detailed imaging and spectral analysis of the available 
ROSAT and ASCA pointed observations are presented. 

The low ionization images are dominated by the bright 
Hii region in the field, not allowing the identification of 
optical emission as emission from the remnant based only 
on morphological arguments. It is the ratio of the [S ii] to 
Ha + [N ii] images which provides some evidence on the 
nature of the emission in this field. We have also detected 
two filamentary structures, unknown up to now, in the 
medium ionization line of [O ill] in the west and the east 
parts of the field observed. The two filaments are found at 
almost opposite sides and display an appreciable degree 
of curvature. The eastern filament is very well correlated 
with the radio emission at 326 and 4850 MHz (Fig. 
while weak radio emission at the location of the western fil- 
ament is only seen in the 1 1 cm maps ( Velusamy & Kundu 
119741 Wendker I1971|l . The strong differences between the 
low and medium ionization lines suggest that significant 
inhomogeneities and density variations in the pre-shock 
medium are present. The Ha, [N ii] and [S ii] lines being 
produced in cooler areas behind the shock front are more 
sensitive to inhomogeneities and density variations than 
the [O III] line which is produced closer to the shock front 
(e.g. Hester il987J . The presence of [Oiii] emission at only 
certain areas of a remnant seems to be a common charac- 
teristic of several remnants like CTB 1 (Fesen et al. I1997|l 
and G 17.4-2.3 (Boumis et al. lM)2| . 

The optical spectrum obtained from area I is a typ- 
ical spectrum from a complete shock structure suggest- 
ing a column density behind the shock front of ~10^^'^ 
cm-2 (e.g. Cox & Raymond HHSl Raymond et al. lTn55|l . 
However, the large [Oiii]/H/3 ratio in area II can only 
be explained by an incomplete shock structure where the 
shock is traveling with a speed greater than 100 km s^^ in 
an interstellar cloud of medium density like in area I. The 
swept-up column density is estimated around 10^^ cm~^ 
(Raymond et al. 11988(1 . The large sulfur line ratios mea- 
sured in the long-slit spectra are indicative of low electron 
densities which in turn imply low pre-shock cloud densi- 



ties (Fesen & Kirshner ll980|l . This is in agreement with the 
cloud density estimates based on the relative line intensi- 
ties of the optical spectra and the modeling of Raymond 
et al. ()1988|l . Currently, the long-slit spectra from area III 
suggest emission from a photoionized nebula. However, the 
complex environment of G 82.2-1-5.3 does not permit an 
unambiguous determination of the nature of the optical 
emission in this area and radio observations at different 
frequencies are required to resolve this issue. 

In an effort to estimate the column density towards 
G 82.2+5.3, based on the optical data, we use the statis- 
tical relation of 

Nh = 5.4 (±0.1) ■ 10^1 E(B - V) cm-2, (1) 

given by Predehl and Schmitt (1995). The estimated neu- 
tral hydrogen column density given the observed values 
of E(B-V) is then in the range of 1-4 x 10^^ cm~^, al- 
though larger column densities cannot be excluded given 
the uncertainties in the optical extinction. Furthermore 
the locations which are accessible for Nh and E(B-V) are 
different. We formally know Nh only for the central part 
of the remnant. Nevertheless, the agreement with the ab- 
sorbing column density derived from the X-ray spectrum 
(4.7x 10^^ cm~^) is not too bad. The full range of X-ray 
absorption of 2.2 - 5.0x10^^ cm~^ corresponds to 0.4 < 
E(B-V) < 0.9. The E(B-V) color excess in the south is 
0.94 and would formally place area III at a larger dis- 
tance than G 82.2+5.3. However, the low counting statis- 
tics from area II and the complex environment do not 
allow to draw secure conclusions on the relative distance 
between the remnant and the emission in the south (area 
HI). Additionally, the upper limit on the X-ray absorp- 
tion does not exclude such a large value of E(B-V) and a 
correspondingly large distance. 

The X-ray emission consists of a patchy disk slightly 
elongated in the south-north direction, dominated by a 
bar-like emission close but clearly off-set from the geo- 
metric center. The remnant is not obviously limb bright- 
ened but a weak shell of emission might be present in the 
eastern hemisphere. Given the available imaging and spec- 
tral X-ray data, we proceed to derive some fundamental 
parameters of G 82.2+5.3 starting with the Sedov rela- 
tion. In this framework the remnant is considered to be 
in the adiabatic phase and estimates for the distance d, 
linear diameter D, age r, luminosity Lx, shock velocity 
Vs, pre-shock matter number density np and the amount 
of swept-up matter Mgu can be derived from the angular 
diameter, the emission measures and the spectral tem- 
peratures (e.g. Pfeffermann et al. 1991). It is interesting 
to note that a reasonable solution of the Sedov relation 
cannot be obtained adopting the temperature of the soft 
component (0.2 keV; see §4.2) as the temperature at the 
primary shock front. However, the temperature of the hard 
component (~0.6 keV), equivalent to a shock velocity of 
~720 km s^^, does provide a viable solution of the Sedov 
relation. Assuming a distance of 1.6 kpc, the explosion en- 
ergy would be 0.17-10^^ erg, the age would be 13500 yr, 
the ISM density would be 0.012 cm~'^, and the swept up 
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Fig. 5. The raw soft X-ray emission detected in the 8 ks ROSAT pointed observation is shown in this figure. The 
brightness scales hnearly from 2 to 24 counts/pixel and the image has been smoothed with a Gaussian filter with (t=8'. 
The radio contours at 326 MHz scale also also linearly from 0.015 to 0.09 Jy/beam (Rengelink et al. I1997|l . 



mass would be 24M0. If, on the other hand, we adopt the 
nominal value of the explosion energy (1-10^^ erg), then 
the age would become 26700 yr, the current radius would 
be 50.4 pc, the ambient density would be 0.009 cm^"^, the 
swept up mass would be 144 Mq , and the remnant should 
be ~3.3 kpc away from us. 

The assumption of pressure equilibrium between the 
interstellar clouds and the hot gas can provide an inde- 
pendent estimate of the pre-shock cloud density. For a 
shock velocity of ~100 km s^^ into the clouds and an am- 
bient density of 0.012 cm^'^, a pre-shock cloud density 



around 1 cm ^ is estimated, which is compatible with the 
estimates based on the long-slit optical spectra (see §3.3). 

With the above data, we can discuss the evolutionary 
status of the remnant in some more detail. Using the cal- 
culations of Cox et al. t l999J and Cioffi et al. (1988) . it is 
found that the cooling times exceed 1.3T0^ yr, i.e. a factor 
of 10 longer than the age estimated above. Furthermore, 
with a radius of 50 arcmin at a distance of 1.6 kpc the 
mean expansion velocity of the remnant over that very 
long age is just 185 km s~^, and with such a low expan- 
sion velocity one might wonder about the X-ray emission 
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Fig. 6. The brightness profile (crosses) of G 82.2+5.3 derived from the ROSAT data is shown in this figure. The open 
squares mark the exponential best fit with some limb brightening, while the open triangles mark the power-law best 
fit. The reader is referred to S4.1 for more details. 



which has a temperature of at least 0.2 keV. Alternatively, 
the remnant might be at a greater distance, say by a fac- 
tor of two to three, which might fit the X-ray temperature 
but it is at the edge of the distance estimates via Nh and 
E(B-V). Furthermore the measurement of overabundant 
Si, which indicates that ejecta matter at fairly high tem- 
peratures still exists, does probably not favour a very old 
remnant. Of course it is not excluded that G 82.2+5.3 has 
reached the radiative phase, despite the very low ambient 
density, which we stress is a direct result of the X-ray mea- 
surements and is independent of any assumption about the 
actual evolutionary status of the remnant. In addition, we 
have examined HI data from the Canadian Galactic Plane 
survey (CGPS - Taylor et al. lMT^ . Unfortunately, due 
to its high galactic latitude of ^^5° only emission from the 
eastern areas of the remnant was recorded in this survey. 
So, at least for this part of the remnant, we were not able 
to identify any features that would point to an HI shell. 

The low ambient density may result from the inter- 
action of the stellar wind of the progenitor star with its 
environment since a typical wind blown bubble is charac- 
terized by a radius of ~30 pc, a density of ~ 0.01 cm~^ and 
a hot interior (e.g. Castor et al. 119751 Weaver et al. ll977|l . 



The low surface X-ray brightness of the annulus around 
G 82.2+5.3, if associated with it, would support the sug- 
gestion of the wind blown bubble in which the progenitor 
star of G 82.2+5.3 exploded and its remnant is expanding. 

The ASCA spectra of G 82.2+5.3 show that the central 
X-ray emission is of thermal origin, while the ROSAT im- 
ages show a azimuthally averaged brightness profile which 
is consistent with a slow exponential decrease of the mat- 
ter density by a factor of about four from the center to the 
edge. In principle, the evaporating cloud model proposed 
by White and Long (1991J or the model of Maciejewski 
et al. H1999|l . where a remnant evolves in a stratified 
medium, can be of relevance to W 63. This model could 
probably explain the patchiness and the cloudy struc- 
ture. Furthermore, heat conduction may apply in order 
to smooth the temperature distribution in the interior of 
the remnant and make it look isothermal (e.g. Chevalier 
1999). However, with the currently available data we can- 
not claim the existence of such mechanisms. 

At last, there is an interesting aspect which comes 
to mind after the analysis of the ROSAT PSPC spec- 
tra. There may be the possibility that the X-ray emitting 
gas we see is just a minor fraction of the total matter in- 
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Fig. 7. The ASCA GIS2/3 data are shown in the left figure together with the best fit vmckal model, while the spectra 
from the SIS detectors are shown in the right figure with the same model (see §4.2). 



side the remnant. If there were gas at temperatures below 
about 0.05 keV, it would have been missed by the X-ray in- 
struments, even if the emission measure would be higher 
by two to three orders in magnitude, basically because 
the intersteller absorption column density is so high, and 
this holds for all the mixed morphology remnants listed 
(Nh > 3.8x10^1 cm-2) by Rho & Petre (,1998.1 . 



6. Conclusions 

We have obtained deep flux calibrated images of 
G 82.2+5.3. The low ionization images are dominated by 
the H II region in the field but the medium ionization line 
of [O III] provides the first clear picture of the remnant. It 
displays an elliptical shape with filamentary structures in 
the east and west areas. The correlation between the op- 
tical and radio emission supports their association. Long- 
slit spectra suggest both complete and incomplete recom- 
bination zones, shock velocities around 100 km s^^, and 
low electron densities. The X-ray emission region appears 
to largely be embedded in the optical and radio emission. 
The X-ray surface brightness is quite patchy, missing ob- 
vious limb brightening and is dominated by a bright bar- 
like emission region which is off-set from the geometric 
center by ~9'. The X-ray emission is definitely thermal 



and requires two temperatures of 0.2 keV and 0.63 keV. 
The bright bar region shows overabundant Mg, Si and Fe, 
which might indicate still radiating ejecta matter. The 
azimuthally averaged radial surface profile is consistent 
with the matter density changing with distance r from 
the center oce"'^/'^" with a characteristic angular length of 
36'; but it is also consistent with an density profile. 

The matter inside the remnant is quite likely structured 
like a porous cloudy medium. The average matter density 
derived from the X-ray surface brightness is ~0.05xdj~g'^ 
and r^0.03xdj^g'^ for the two spectral components, respec- 
tively. It could be slightly higher with filling factors sig- 
nificantly less than 1. From the X-ray and optical mea- 
surements we have determined some parameters of the 
SNR, but the evolution of the remnant into the mixed- 
morphology state is still difficult to understand. If we ap- 
ply the 'Sedov analysis' approach, which would assume 
the remnant to be in the adiabatic phase, the remnant 
is at a distance between 1.6 kpc and 3.3 kpc with an 
age of 13.5 and 26.7 kyrs. These results appear reason- 
able, particularly because the distance is consistent with 
that determined by Rosado and Gonzales ((1981) . But the 
remnant shows neither the expected density profile nor a 
strong limb brightening for a constant ambient density on 
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Fig. 8. Best fit and error contours (Icr to 4(t in steps of la) for the temperature (kT) and the interstellar absorption 
column density Nh for one-temperature vmekal models to the ROSAT PSPC spectrum (lower left) and the ASCA 
GIS/SIS spectra (upper right). Note that there is essentially no overlap, indicating the need for more advanced models. 



which the 'Sedov analysis' is based. In addition, we de- 
termine the temperature only at the center and not at 
the limb which is what is required for the application of 
the 'Sedov analysis'. On the other hand, the assumption 
that the remnant is in the radiative phase entering the 
shell-forming phase is in conflict with the very low den- 
sity of the X-ray emitting plasma and the excessive long 
radiative cooling time involved. Although the low mat- 
ter density itself does not completely rule out that the 
remnant has reached the radiative phase it would imply 
an even lower density, greater age and much larger dis- 
tance at the edge of the upper limits obtained from Nh 
and E(B-V). These conclusions arc primarily based on the 
low density of the X-ray plasma but the optical filaments 
demonstrate the existence of matter of pre-shock densities 
of one or a few atoms cm"'^ and shock velocities of 100 
km s~^. Whether these filaments are located around the 
perimeter of the remnant or inside of it is not clear. But 
matter at these densities and velocities of >100 km s~^ 



should emit soft X-rays or EUV radiation. The fact that 
we do not see such radiation puts an upper limit on the 
temperature and the emission measure, or electron den- 
sity, of this component in order to stay 'undetected' in 
the PSPC observations which depend on Nh- Adopting 
the observed value of Nh the pre-shock electron density 
could be as high as 1 cm~'^ for log(T) = 5.8 or kT = 55 eV. 
Interestingly, the radiative cooling time is 25 kyrs, and 
an HI shell might not have formed yet. We have not found 
any evidence for an HI shell, but this is not conclusive 
as the full area around the remnant has not been imaged 
so far. Even higher electron densities up to 20 cm~^ with 
log(T) = 5.5 are possible for the soft X-rays/EUV to es- 
cape detection but the radiative cooling time gets as short 
as 10 kyrs, and the probability to have caught a remnant 
in such a brief transition phase gets increasingly lower. 
If most of the dense matter were predominantly located 
around the perimeter of the remnant it might represent 
the shock heated matter of parts of the inner edge of a stel- 
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lar wind-blown bubble, which might still be larger than 
the remnant, indicated by the low- level X-ray emission 
surrounding the remnant as seen in the ROSAT All-sky 
survey. If the supernova had a stellar wind-blown bubble, 
the low electron density derived from the X-ray measure- 
ments is easy to understand as it matches the prediction of 
densities in such bubbles. The temperature and the X-ray 
surface brightness distribution observed in the remnant 
are then a consequence of the temperature and matter 
density distributions in the bubble prior to the supernova 
shock, both of which change substantially with distance 
from the progenitor star (Weaver et al. 11977(1 . 
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